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Abstract 
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Previous  experiments  on  plasma-propellant  interaction  have  shown  enhanced  bum  rates.  The 
temperature  of  the  propellant’s  bed  may  play  a  role  in  the  bum  rates  under  plasma  injection.  In  this  report, 
a  set  of  experiments  were  first  conducted  on  the  plasma-flow-field  to  determine  the  spatial  and  temporal 
distribution  of  the  plasma  temperature,  pressure,  number  density,  and  velocity.  The  experiments  revealed  a 
decreasing  plasma  pressure,  plasma  temperature  and  plasma  number  density  as  plasma  is  leaving  the 
capillary  source  and  expands  in  air.  The  plasma  jet  velocity  2  inches  from  the  source  exit  was  found  to  be 
about  1300  m/s.  Following  characterization  of  the  plasma-flow-field,  a  set  of  experiments  were  conducted 
on  JA-2  solid  propellant  with  controlled  bed  temperature.  Experiments,  with  proposed  models,  revealed  a 
functional  form  that  includes  the  bed  temperature.  Increased  bum  rates  were  observed  with  increased  bed 
temperature,  which  appears  to  follow  a  power  law.  A  model  for  the  bum  rate  is  proposed,  which  includes 
the  bed  temperature  and  has  the  BR  =  A  Pa  (  T/  T  ambient  )  b.  Further  investigation  on  plasma 
parameters  helped  in  modifying  the  model  to  include  plasma  radiative  heat  flux  and  plasma 
number  density  in  the  proposed  model. 

I.  Introduction 


Electrothermal  plasma  has  been  studied  as  a  viable  replacement  for  solid  propellant  ignition. 
The  studies  were  not  only  conducted  for  possible  ignition  using  plasmas,  but  also  as  means  of 
controlling  and  enhancing  bum  rates  of  propellants.  With  plasma  augmenting  the  propellant’s  bum 
rate,  a  mix  between  plasma  electrical  energy  and  chemical  energy  of  the  propellant  may  provide  a 
possible  means  of  achieving  high  enthalpy  flows.  Such  high-enthalpy  flows  may  provide  better 
interior  ballistics  and  higher  muzzle  velocities  for  hyper-velocity  launchers  [1-8].  Previous  work 
has  shown  that  the  plasma  ignited  propellant  bum  increases  surface  erosion  compared  to  that  of 
conventional  ignition  [1-5].  Those  experiments  studied  the  effects  of  plasma  pressure  impact  on 
the  bum  rate  of  propellant  when  the  propellant  was  at  ambient  conditions.  No  propellant  bed 
temperature  sensitivity  has  been  examined,  as  of  yet’  when  the  propellant  is  ignited  by 
electrothermal  plasmas. 

This  work  will  determine  if  any  propellant  temperature  dependence  exists  in  the  bum  rate  and, 
if  it  exists,  the  functional  form  of  that  dependence.  A  series  of  experiments  was  performed  on 
NCSU’s  PIPE  [2,4,5,  9,10]  ETC  device  that  impacted  plasma  onto  a  heated  propellant  bed  (JA2). 
Experiments  were  conducted  using  an  incomplete  “extinguished”  bum  technique  to  provide  an 
evaluation  of  the  erosive  bum.  These  plasma  impacts  resulted  in  propellant  mass  loss  that  could 
be  correlated  to  bum  rate.  If  propellant  bed  temperature  significantly  impacts  the  bum  rate  of  the 
propellant  then  this  would  be  shown  in  the  experimental  data  returned  from  the  mass  loss  as 
compared  to  bum  rates  at  ambient  temperature.  The  experiments  also  required  a  pressure 
difference  so  the  propellant  was  placed  at  1,  2  and  3  inches  from  the  barrel.  This  allowed  for  a 
correlation  based  on  pressure  and  propellant  bed  temperature.  Additionally,  the  radiative  heating 
effect  of  the  plasma  on  the  propellant  bed  could  be  evaluated  from  the  plasma  temperature  change 
farther  away  from  the  barrel. 

The  main  goal  of  the  work  is  to  determine  a  functional  equation  that  can  be  used  with 
common  diagnostics  to  determine  the  impact  of  the  bed’s  temperature  on  bum  rates.  Pressure  is 
commonly  measured  in  expansion  chambers  and  is  relatively  easy  to  obtain.  As  a  result,  the 
functional  form  will  have  pressure  and  bed  temperature  dependence  and  the  bum  rate  with  these 
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dependencies  can  be  calculated.  The  final  goal  of  all  current  propellant  studies  is  the  development 
of  a  plasma-sensitive  propellant  that  will  only  combust  when  burned  with  plasma.  This  will  be  a 
safer  propellant  that  will  not  combust  under  fire  or  by  any  other  accidental  means.  This  work 
continues  the  efforts  to  develop  such  a  plasma  sensitive  propellant  for  the  viable  use  of  ETC 
launchers. 


Many  researchers  have  studied,  in  detail,  the  interaction  of  plasma  with  propellants. 
Electrothermal  plasmas  are  used  because  they  are  characterized  by  their  high-density  and 
relatively  low-temperature  (1-1  OeV)  such  that  ET  plasma  resembles  a  blackbody  [7- 16]. 

The  majority  of  the  research  at  the  ARL  has  been  performed  in  order  to  determine  the  effects 
of  plasma  temperature,  pressure,  and  radiation  effect  [16-21],  In  general,  these  are  compared 
with  the  conventional  black  powder  ignition  means  in  a  large  assembly  that  simulates  how  a  full 
propellant  packed  charge  would  bum.  These  provide  a  good  indication  of  how  a  real  cannon 
device  would  be  operated  and  an  idea  of  the  geometry  of  the  bum  propagating  through. 
However,  these  closed  bomb  devices  don’t  lend  themselves  well  to  the  study  of  optimization 
because  of  the  burning  and  pressure  due  to  other  grains.  Individual  grain  study  is  useful  for 
determining  the  actual  effect  of  plasma  on  the  propellant  and  gives  an  idea  as  to  the  improvements 
in  packing  geometry  or  even  chemical  composition  to  give  the  desired  bum  characteristics. 
NCSU’s  plasma  engineering  group  has  ETC  devices  that  are  ideally  suited  for  single  propellant 
grain  plasma  experiments  and  as  a  result  are  contracted  by  ARL  to  perform  propellant  studies. 
Previous  work  has  shown  geometry  effects  by  rotating  the  propellant  to  change  the  plasma  impact 
angle.  The  study  revealed  that  normal  impact  results  in  the  highest  bum  rate. 

The  traditional  form  of  the  combustion  equation  comes  from  chemical  kinetics.  For 
example  a  bimolecular  reaction  of  the  type 

A  +  B  —>  M  +  N 


The  rate  of  reaction  is  proportional  to  the  probability  of  collision  in  the  molecules.  This 
probability  is  determined  by  the  number  of  moles  per  unit  volume  or  the  concentrations.  Hence 


dNjL=  d(VCA)_  dNa 
dt  dt  dt 


djy£s) 

dt 


=  kVCACB 


where  V  is  the  volume  and  k  is  the  reaction  rate  constant. 


Performing  differentiation  and  dividing  through  by  a  constant  volume  yields 


dC ± 
dt 


d£ \ 
dt 


=  kCACB 


If  the  reaction  proceeds  as  the  stoichiometric  equation  states 


aA  +  bB  — »  mM  +  nN 
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then  the  rate  of  change  of  the  concentration  is  related  to  those  remaining  and  can  be  written  in  the 
form  [22] 


d£, 

dt 


=  knC 


n 

A 


The  concentration  can  be  written  as  density  times  the  mole  fraction  and  taking  the  full 
derivative  of  the  concentration  yields 


Dividing  by  a  constant  density  and  solving  for  the  change  in  mole  fraction  gives 


ds 

—  =  -kpn-'en 
dt 


The  mole  fraction  lost  due  to  combustion  is  proportional  to  the  bum  depth  of  the  material  and 
density  times  the  mole  fraction  is  proportional  to  pressure.  This  allows  the  equation  to  be  written 
in  the  conventional  form  of  the  Saint  Roberts  Equation  [5]: 


BR  =  bPa 


with  P  in  psi  where  a  takes  care  of  the  combination  of  n  and  n-1 .  The  value  for  the  exponent  n 
was  determined  experimentally  for  the  modified  double  base  propellant  JA2  to  be  n=1.3755.  The 
b  was  determined  to  be  8.865x1  O'5.  This  value  is  higher  for  plasma  impacted  propellant  than  for 
the  conventional  black  powder  combustion  n=0.889  and  b=203.2xl0'5  [4],  This  results  in  higher 
bum  rates  for  an  equivalent  pressure  (because  the  value  of  b  is  constant  for  each  propellant) 
which  would  change  the  ballistics  of  a  projectile  launched  by  the  propellant  (and  is  necessary  to 
know  in  the  case  of  cannons).  Bourham  also  showed  experimentally  that  the  bum  rate  is  not 
dependent  solely  on  pressure  but  is  also  coupled  with  the  radiation  heat  flux  from  the  plasma 
temperature.  This  was  done  by  comparing  effects  of  single  and  double  exposed  samples  [4].  The 
samples  that  were  exposed  twice,  once  on  each  face  of  the  propellant,  showed  a  higher 
gasification  rate.  The  conclusion  was  drawn  that  this  gasification  rate  was  due  to  the  in-depth 
reactions  caused  by  the  radiation  flux  of  the  plasma. 

It  has  also  been  shown  that  the  only  augmentation  the  plasma  gives  to  the  propellant  bum  rate 
is  during  the  actual  plasma  impact  and  there  are  no  secondary  bum  rate  augmentations.  This  was 
shown  by  Birk  [8]  who  also  did  some  experiments  with  chilled  propellants.  His  results  with  the 
M30  propellant  showed  that  the  bum  rate  at  -22  C  was  significantly  lower  than  those  samples 
exposed  at  22  C.  Even  though  he  didn’t  have  any  comparative  data  for  the  JA2  propellant,  one 
can  conclude  that  a  similar  effect  should  take  place  with  the  JA2  propellant. 

If  this  temperature  effect  had  a  considerable  impact  on  the  bum  rate  of  a  propellant,  then  this 

effect  must  be  quantified.  The  artillery  officer  in  charge  of  the  weapon  and  fire  guidance  must 
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know  how  the  environment  is  going  to  impact  the  range  and  the  exterior  ballistics.  If  the  weapon 
is  deployed  into  an  environment  similar  to  Siberia  or  to  an  environment  like  the  Sahara  Desert  the 
artillery  officer  may  not  have  time  to  calibrate  the  devices  using  a  trial  and  error  method. 

Information  about  the  firing  of  the  propellant  should  be  available  in  order  to  determine  if  the 
environment  will  affect  the  firing  and  to  what  extent  the  environment  will  impact  the  propellant. 

For  this  reason,  it  is  proposed  to  modify  the  bum  rate  equation  to  include  propellant  bed 
temperature  dependence  in  the  form 

BR{P,T)  =  aPbf(T) 

All  of  the  previous  experiments  to  determine  the  constant  a  and  exponent  b  were  done  at  an 
ambient  temperature  assumed  to  be  20  °C,  but  the  propellant  temperature  was  not  monitored. 
This  work  will  determine  if  the  propellant  temperature  does  have  a  significant  impact  on  the  bum 
rate  of  plasma  ignited  propellant  and  if  so,  what  the  pressure-temperature  dependent  bum  rate 
would  be.  It  is  assumed  that  the  bum  rate  dependence  on  pressure  will  remain  as  a  power  law  but 
it  is  possible  that  the  temperature  dependence  could  follow  a  power  or  an  exponential  form 


BR  =  aPbTc 
or 

BR  =  aPb  exp (cT) 


The  bed  temperature  may  also  be  expressed  as  a  ratio  between  the  raised  temperature  and  the 
temperature  at  ambient: 


BR 


aPb 


f  y  ' 
Jo  J 


where  T0  is  the  ambient  temperature. 

The  expectation  is  that  there  will  be  propellant  temperature  sensitivity  in  the  bum  rate.  Beyer 
[7]  already  showed  some  effect  of  cooling  of  M30  propellant  and  this  side  note  found  in  his  work 
must  be  expanded  on.  However,  it  makes  sense  logically  that  the  temperature  of  the  propellant 
bed  would  have  an  effect  on  bum  rate  because  as  the  temperature  changes  in  the  bed,  so  does  the 
energy  of  the  bed.  Molecules  that  are  excited  are  more  likely  to  dissociate  from  each  other  when 
ignited,  thus  increasing  the  bum  rate. 


II.  Plasma  Flow-Field  Characterization 


Plasma  flow-field  characterization  and  distribution  of  plasma  parameters,  both  spatial  and 
temporal,  gives  an  accurate  picture  on  the  environment  where  the  propellant  bed  is  placed  within 
during  exposure  to  electrothermal  plasma.  With  controlled  bed’s  temperature  to  evaluate  the 
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sensitivity  of  the  propellant  to  the  plasma  jet  at  elevated  temperatures,  it  is  necessary  to  determine 
the  plasma  parameters  in  the  plasma- flow  field.  Plasma  pressure  and  temperature  of  the  plasma 
jet  outside  the  electrothermal  source  must  be  known  by  direct  measurements.  These 
measurements  can  help  determining  the  momentum  and  energy  terms  at  any  chosen  location,  and 
thus  the  plasma  impact  on  the  propellant  bed  can  be  expressed  in  terms  of  the  plasma  parameters 
at  chosen  location.  Before  conducting  temperature  sensitivity  studies,  the  plasma-flow-field 


experiments  were  conducted  and  the  pressure,  temperature,  density  and  velocity  distribution  in  the 
flow-field  were  measured.  A  compact  multi-sensor  probe  has  been  designed  and  constructed  to 
measure  spatial  and  temporal  plasma  parameters.  These  parameters  are  the  total  and  Static  plasma 
pressure,  plasma  temperature,  plasma  number  density  and  plasma  Velocity.  All  measurements 
were  conducted  at  atmospheric  pressure  by  fusing  the  ET  source  with  aluminum  fuses.  The  ET 
plasma  is  thus  jetting  into  open  air  inside  of  a  cubical  chamber.  The  plasma  jet  expands  out  of  the 
ET  source  with  a  measured  expansion  angle  of  approximately  23°  from  the  ET  source’s  axis. 

1M  The  Compact  Multi-Sensor  Probe: 

The  compact  multi-sensor  probe  [23],  as  shown  in  the  Fig.  1  below,  has  two  absolute  pressure 
transducers  (Kistler  Model  60 IB)  and  three  fiber  optic  cables.  The  sensor  head  is  arranged  such 
that  the  two  pressure  transducers  are  close  to  each  other  with  one  pointing  towards  incoming 
plasma  flow  to  measure  total  pressure,  and  one  situated  upright  to  measure  the  static  pressure. 


Compact  probe  head: 
2  pressure  transducers 
and  3  fiberoptics 


Fig.  1  Schematic  of  the  concept  of  multi-sensor  probe  equipped  with  two  pressure  transducers 
and  three  fiber  optic  cables.  Positioning  of  the  pressure  transducers  is  shown  at  right. 

To  calibrate  the  pressure  transducers,  a  shock  tube  was  designed,  constructed  and  used  for 
calibration.  It  consists  of  a  high-pressure  chamber  and  a  barrel  section,  separated  by  a  burst  disc. 
The  high-pressure  chamber  is  connected  to  a  compressed  nitrogen  cylinder  and  set  to  the  desired 
pressure.  The  valve  is  shut  and  the  device  is  triggered  abruptly  via  a  plunger,  which  in  turn  forces 
the  burst  disk  to  burst  and  thus  forwarding  a  shock  wave  down  the  test  barrel.  Pressure  signals 
are  recorder  on  a  LeCroy  9310  digital  storage  oscilloscope.  It  was  determined  that  the  relation 
between  the  true  pressure  versus  measured  pressure  is  linear  and  is  expressed  by  Ptrue  =  1.03 
P measured  -  0.46.  The  three  fiber  optic  cables 
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are  arranged  as  a  bundle,  with  one  fiber  optic  interfaced  to  an  optical  multichannel  analyzer  for 
time-integrated  optical  emission  spectra  measurements,  and  two  other  fiber-optics  each  is 
interfaced  to  a  monochromator  and  a  photomultiplier  for  time-resolved  spectral  measurements. 
Monochromators  were  adjusted  to  521.82  nm  and  570.02  run  wavelengths  (copper  lines).  Optical 
emission  spectroscopy  [24-26]  was  used  to  measure  the  plasma  temperature,  assuming  plasma  is 
at  LTE,  and  plasma-density  from  optical  lines  broadening.  Three  axial  locations  were  chosen  for 
measurements,  2.4,  3.4  and  4.4  inch  from  source  exit.  At  each  axial  location,  three  radial  locations 
are  defined  to  measure  distributions  within  the  diverging  plasma  jet.  Figure  2  bellow  shows  the 
source  exit  and  the  locations  of  compact  multi-sensor  probe  measuring  positions. 


Fig.  2  Schematic  showing  selected  locations  to  obtain '  spatial  distribution  of  electrothermal 
plasma  parameters  down  stream  from  the  source  exit  and  within  a  20°  angle. 

Selection  of  the  measuring  locations  is  based  on  the  plasma  jet  divergence  as  the  plasma  leaves 
the  source.  The  angle  was  measured  by  placing  a  witness  plate  in  the  plasma  stream  along  the  axis 
of  the  device.  The  nine  selected  location  covers  the  solid  angle  of  divergence  and  provides  a 
detailed  information  on  the  spatial  distribution  of  the  plasma  jet  parameters. 

II.2  Plasma-Flow-Field  Distribution: 

Plasma  parameters  at  2.4,  3.4  and  4.4  inches  from  the  source  exit  were  measured,  and  at  three 
radial  locations  for  each  axial  position.  Figures  3  and  4  show  the  pressure  distribution  from  the 
source  exit.  Figure  3  shows  the  results  of  the  peak  total  pressure,  where  it  is  obvious  that  the  axial 
profile  is  approximately  flat  for  all  radial  locations,  while  the  radial  profiles  are  of  a  top-hat 
distribution.  The  measurements  imply  that  for  axial  locations  between  2  and  4  inches,  the  plasma 
jet  undergoes  an  essentially  isentropic  expansion,  while  the  edge  of  the  jet  is  subjected  to  mixing 
process  with  the  surrounding  gas. 

Figure  4  shows  the  distribution  of  the  peak  static  pressure,  where  the  general  trend  indicates  a 
decrease  in  the  static  pressure  as  the  jet  expands.  Of  course  the  largest  static  pressure  is  at  the 
location  closest  to  the  source  exit. 
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Fig.  3  Peak  total  pressure  in  the  expanding  plasma  jet 


Spatial  distribution  of  the  peak  static  pressure  (psig)  of  the  plasma  jet  as  measured  from  the  source  exit 


Fig.  4  Peak  static  pressure  in  the  expanding  plasma  jet 


Of  importance  is  the  distribution  of  the  plasma  temperature  and  density.  The  spatial 
distribution  of  the  plasma  temperature  is  shown  in  Fig.  5.  These  temperatures  were  obtained  from 
time-integrated  optical  emission  spectroscopy  using  copper  lines.  The  plasma  temperature  drops 
quickly  from  2.1  eV  right  at  the  source  exit  (axial  distance  =  0)  to  approximately  0.5  eV  at  2 
inches  down  stream.  Further  into  the  stream,  one  may  observe  that  the  plasma  jet  is  approximately 
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isothermal.  With  additional  data  from  time-resolved  spectral  measurements  of  copper  lines,  it  was 
obvious  that  the  plasma  temperature  inside  the  jet  is  0.4  +/-  0.16  eV.  A  slight  increase  in  the 
plasma  temperature  was  observed  at  4.4  inches  from  the  source  exit.  Such  temperature  increase  is 
attributed  to  the  expansion  of  the  plasma  jet. 


Fig.  5  Average  plasma  temperature  in’the  expanding  jet 


The  plasma  density  is  another  important  parameter.  It  was  calculated  from  the  peak  static 
pressure  and  the  time-integrated  plasma  temperature.  The  plasma  density  at  spatial  locations  of 
measured  pressure  and  temperature  is  shown  in  Fig.  6.  The  density  decreases  with  increased  axial 
length  as  well  as  increased  radial  position.  This  is  expected  for  an  expanding  plasma  jet  from  an 
electrothermal  source.  It  is  obvious  that  the  plasma  density  follows  the  static  pressure  distribution, 
especially  that  the  temperature  is  almost  same  inside  the  jet  and  thus  the  density  profile  follows 
the  static  pressure  profile. 

The  average  speed  of  the  plasma  jet  was  obtained  from  the  plasma  arrival  time  to  the  sensors 
at  the  chosen  measuring  locations.  The  speed  of  the  jet  is  approximately  constant  and  close  to 
1300  m/s  on  axis.  The  slowing  of  the  plasma  near  the  jet  boundary  is  attributed  to  the  mixing 
processes  with  the  surrounding  air.  For  a  0.5  eV  ideal  gas,  the  sound  speed  is  3100  m/s,  and  thus 
the  measured  speed  indicates  that  the  jet  flow  is  sub-sonic  with  a  Mach  number  of  about  0.4.  The 
time-resolved  spectral  measurements  were  also  used  to  calculate  the  plasma  precursor,  from 
which  the  precursor  plasma  velocity  is  found  to  be  about  3000  m/s.  Fig.  7  shows  the  average 
plasma  speed  in  the  expanding  jet. 
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The  obtained  parameters  are  important  in  deciding  where  to  place  the  propellant  for  a  given 
experiment,  and  what  to  be  expected  in  terms  of  plasma  parameters,  at  any  selected  location. 


Fig.  6  Plasma  density  distribution  in  the  expanding  jet 


Fig.  7  Average  plasma  speed  in  the  expanding  jet 
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HI.  Temperature  Sensitivity  Studies 

Temperature  sensitivity  of  the  propellant  bed  may  alter  the  bum  rates  under  plasma  injection. 
When  the  bed’s  temperature  is  elevated,  it  is  expected  that  the  bum  rates  would  increase  as  a 
result  of  increased  in-depth  heating  of  the  bed.  Also,  the  plasticizer  may  migrate  causing 
alterations  in  the  bum  rates.  A  temperature-control  unit  is  designed  such  that  the  bed’s 
temperature  prior  to  plasma  injection  is  maintained  constant.  The  bed  temperature  can  be  elevated 
up  to  100°C  ,  however,  the  upper  limit  was  set  to  50°C  [27]. 

II1.1  Experimental  Arrangement: 

A  circulating  cooling  unit  was  designed,  constructed  and  tested  for  operation.  Preliminary 
experiments  on  the  circulating  cooling  control  unit  have  shown  the  necessity  of  long  time 
circulation  of  the  fluid  in  order  to  maintain  temperature  stability.  A  modified  design  using  heating 
elements  allowed  for  a  faster  temperature  setting  prior  to  plasma  injection.  Figure  8  shows  the 
experimental  arrangement  using  heating  elements  to  fulfill  the  temperature  sensitivity  studies  with 
elevated  temperatures  between  20  to  50°C. 


Propellant 


Heating  element 


ET  source 


Heating 
I  element  leads 


/I _ 11111118 

mm — 


Adjustable 
sample  holder 


Fig.  8  Schematic  of  the  experimental  setup  using  heating  a  heating  element  to  elevate  the  bed’s 
temperature  prior  to  plasma  injection  on  the  propellant  sample 


The  sample  holder  is  tube-shaped  and  has  a  seat  for  the  propellant  sample,  which  resides  on 
the  heating  element.  The  temperature  of  the  back  surface  of  the  propellant  is  monitored  via  a 
thermocouple.  The  sample  is  thin,  about  l-2mm  and  the  sample  temperature  is  kept  stable  to  the 
desired  temperature  via  the  heating  element  controller.  The  holder  has  a  fiberoptic  cable  to  view 
the  plasma  at  the  plasma-propellant  interface.  The  entire  holder  is  adjustable  and  can  be  placed  at 
any  location  on  the  axis  of  the  plasma  source.  Experiments  are  conducted  at  a  base  pressure  of  20 
Torr  for  extinguished  bum  (interrupted  bum)  such  that  the  bum  is  only  during  the  plasma  pulse 
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length  (100  ps  discharge).  The  initial  and  final  weights  of  the  propellant  sample  are  recorded,  as 
well  as  the  initial  and  final  weights  of  the  plasma  source  liner  and  the  source  electrode.  For  each 
shot,  the  discharge  current  is  fed  to  our  TITAN  2-D  code  [28]  to  calculate  the  plasma  exit 
parameters.  Plasma  temperature  at  the  plasma-propellant  interface  is  calculated  from  optical 
emission  spectroscopy  measurements  using  copper  lines. 

A  block  diagram  of  the  experimental  arrangement  is  shown  in  Fig.  9,  where  PIPE  represents 
the  ETC  experimental  facility  powered  by  a  360pF,  lOkV  high  energy-density  capacitor.  The 
block  diagram  shows  experimental  components  and  data  acquisition  system,  thermocouples  and 
fiber  optics  to  optical  multichannel  analyzer. 


The  sample  holder  was  designed,  essentially,  to  be  a  rectangular  block  with  supporting  rods 
that  could  be  bolted  to  the  bottom  of  the  ETC  chamber  of  the  PIPE  facility.  The  heater  used  was 
a  3/8"  diameter  OMEGA  cylindrical  cartridge  heater  that  was  placed  on  a  shelf  and  held  in  place 
by  an  aluminum  plate.  A  small  copper  tube  was  attached  to  the  holder  by  epoxy  in  order  to  keep 
the  propellant  sample  from  sliding  against  the  pressure  and  provides  a  way  to  maintain  thermal 
contact.  The  initial  temperature  of  the  propellant  bed  was  monitored  by  a  thermocouple  pressed 
against  the  exposed  face  facing  the  plasma.  The  original  method  for  controlling  the  heater  output 
was  an  Omega  Process  Controller  that  controls  the  heater  temperature  using  an  output  relay. 
This  controller  was  chosen  because  it  would  keep  the  heater  temperature  essentially  constant,  and 
could  be  monitored  and  controlled  through  varying  the  voltage  across  the  heater.  However,  two 
units  were  destroyed  in  a  relatively  short  time  because  of  the  high-voltage  high-current  discharge. 
A  large  potentiometer  was  added  to  the  system  and  used  to  power  the  heater.  The  potentiometer 
regulates  the  voltage  on  the  heater  at  a  given  time,  thus  controlling  the  current  and  temperature  of 
the  heater.  This  setup  required  careful  monitoring  of  the  rate  of  temperature  increase  so  that  the 
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experiment  would  be  performed  at  the  time  corresponding  to  the  desired  propellant  bed 
temperature.  It  was  possible  to  obtain  a  pseudo-steady  state  temperature  with  this  device  after  an 
hour  and  a  half  of  operation  because  the  heat  losses  to  the  environment  would  eventually  balance 
out  the  input  heat  flux  to  the  heater. 

III.2  Bed  Temperature  Sensitivity  Results 

Results  of  the  propellant  bed’s  temperature  sensitivity  have  shown  an  increasing  bum  rates 
with  increased  bed  temperature.  For  the  variation  of  bed  temperature  between  23  to  50°C,  the 
best  fit  is  showing  a  power  law  of  the  form:  Bum  rate  BR  =  A,  Tm,  where  the  bum  rate  is 
measured  by  dividing  the  interrupted  erosive  thickness  by  the  plasma  impact  time  (FWHM).  Here 
in  this  equation,  T  is  the  bed  temperature  in  °C,  Ai  is  proportionality  constant  and  m  is  the  power 
exponent.  This  equation  can  be  introduced  to  the  bum  rate  equation  such  that  the  bum  rate, 
including  bed  temperature,  may  be  expressed  by:  BR  =  A  Pa  Tb,  where  the  first  term  (Pa) 
represents  the  conventional  equation  for  the  bum  rate  without  introducing  the  bed’s  temperature. 
A  more  general  equation  may  be  introduced  as:  BR  =  A  Pa  (  T/  Tambient)b,  and  thus  for  T  =  Tambient , 
the  bum  rate  falls  to  the  known  equation.  Bearing  in  mind  that  the  pressure  exponent  ‘a’  with 
plasma  ignition  is  much  higher  than  its  value  for  conventional  ignition.  The  following  graph  shows 
the  experimental  results  for  elevated  temperature  regime  with  the  propellant  situated  2  inches 
from  the  source  exit.  When  the  propellant  is  closer  to  the  source  exit,  it  is  expected  that  the  bum 
rate  will  be  much  higher  because  of  the  pressure  effect  (momentum  effect)  on  the  bum  rate. 

When  expressing  the  bum  rate  by  the  equation  Bum  Rate  BR  =  A  P“,  it  has  been  shown  for 
JA-2  that  the  bum  rate  in  conventional  ignition  is  given  by  BR  =  203.2x1  O'5  P0  889,  while  for 
plasma  ignition  it  was  shown  to  have  a  higher  exponent:  BR  =  8.865x1  O'5  PU76.  As  stated  above, 
to  include  the  effect  of  the  bed  temperature,  it  is  necessary  to  modify  the  bum  rate  equation  to 
include  the  bed  temperature 

BR  (P,T)  =  a  PbXT)  =  a  Pb  Tc,  where  T  is  the  bed  temperature. 

One  may  express  the  bum  rate  in  terms  of  normalized  bed  temperature  in  the  form 


BR  =  a  Pb 


KTo) 


Where  To  represents  the  ambient  (room)  temperature.  This  model  is  more  convenient  to  relate  the 
bum  rate  to  the  ratio  between  elevated  temperature  to  ambient  temperature,  or  could  be 
expressed  as  the  bum  rate  change  for  a  percentage  change  in  the  bed  temperature.  The  advantage 
of  this  model  is  that  it  directly  reduces  to  the  Saint  Roberts  bum  rate  at  ambient  temperature. 
Figure  10  shows  experimental  results  obtained  for  JA-2  bum  rates  versus  increased  bed 
temperature  for  one  data  set  at  2-inches  from  the  electrothermal  source  exit.  It  is  obvious  that  the 
data  easily  fits  a  power  law  with  a  2.4767  exponent. 
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Fig.  10  Bum  rate  of  JA-2  with  increased  bed  temperature.  Sample  placed  2  inches  from  the  exit  of 
the  ET  source.  Data  fit  to  power  law  with  2.4767  temperature  exponent 

The  bum  rate  data  for  all  selected  axial  locations  of  1 ,  2  and  2.75  inches  are  shown  in  Fig.  1 1 
as  a  function  of  the  normalized  temperature  (T/T0).  The  bum  rate  data  shows  that  the  1-inch  shots 
have  the  highest  bum  rate  and  the  2-inch  the  lowest.  Initially  one  may  think  that  the  2.75-inch 
data  should  be  the  lowest  and  these  data  are  in  error.  However,  referring  to  the  pressure  data 
obtained  by  the  compact  multi-head  probe  as  previously  shown  in  section  II,  it  is  shown  that  total 
plasma  pressure  actually  increases  slightly  between  2  and  3  inches  as  a  result  of  momentum 
transfer  and  expansion. 

Since  the  pressure  of  these  data  was  not  directly  measured,  a  method  of  evaluating  the 
pressure  is  to  compare  the  source  exit  pressure  with  previously  measured  pressures.  The  plasma 
diagnostic  code  TITAN  was  run  to  determine  the  barrel  exit  pressure  at  atmospheric  conditions 
and  then  this  value  was  fitted  to  the  values  obtained  by  compact  probe  measurements  (as 
previously  described  in  section  II).  TITAN  is  a  pseudo  2-D  code  that  calculates  the  plasma 
parameters  at  the  source  exit  both  spatially  and  as  a  function  of  time.  TITAN  was  also  used  with 
base  pressure  of  20  Torr  (the  experimental  conditions)  and  compared  to  compact  probe  results  to 
see  if  the  chamber  pressure  had  any  significant  effect  on  plasma  pressure.  TITAN  revealed  that 
there  is  no  significant  change  between  the  plasma  expanding  at  20  Torr  or  expanding  at  760  Torr 
[6],  The  plasma  barrel  pressure  was  approximately  1 1,400  psi  independent  of  the  initial  chamber 
pressure.  This  value  was  fitted  to  compact  probe  data  at  2,3,  and  4  inches  to  approximate  the 
plasma  pressure  at  corresponding  distances.  The  pressure  values  are  shown  in  Table  1  below. 


Table  1:  Plasma  Pressure  (Total) 


Distance  (inches) 

- ^ - 6— _ _ _ _ 

Pressure  (psia) 

1 

1403.6 

2 

295 

3 

304 

14 


Fig.  1 1  Bum  rate  of  JA-2  with  increased  bed  temperature  for  all  samples’  positions.  Normalized 
temperature  (T/Tc)  is  used  in  this  graph 


Matlab’s  nonlinear  least  squares  function,  lsqnonlin,  was  used  to  determine  the  unknowns  in 
the  modified  bum  rate  equation  with  values  shown  in  Table  2. 
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Thus,  the  modified  bum  rate  equation  may  be  written  as 

BR  =  0.355P0'528 


f  T  V-36 


\Toj 


JA  2 


This  equation  fits  well  the  data  set  obtained  at  all  measured  locations,  as  seen  in  Fig.  12. 


Fig.  12  Bum  rate  of  JA-2  with  increased  bed  temperature  for  all  samples’  positions.  Normalized 
temperature  (T/Tc)  is  used  in  this  graph.  Data  fitting  used  normalized  bum  rate  equation. 

However,  it  can  be  seen  that  the  2-inch  data  is  not  well  fitted.  This  indicates  there  is  another 
factor  that  is  not  taken  into  account  in  the  above  equation.  The  equation  should  reduce  to  the 
numbers  obtained  by  Zaghloul  et  al  [5]  when  the  propellant  bed  reaches  ambient  temperature. 
This  was  not  the  case.  The  bum  rates  obtained  in  this  experiment  at  propellant  ambient 
temperatures  were  much  higher  than  those  obtained  previously  at  similar  pressures.  The  difference 
is  attributed  to  the  way  the  nonlinear  solver  fits  the  function.  A  quick  calculation  at  1 1,400  psi 
and  propellant  bed  at  ambient  still  shows  an  increased  bum  rate  with  plasma  injection  over 
conventional  means.  With  plasma  injection  the  bum  rate  at  ambient  temperature  is  49.2  cm/sec 
verses  8.2  cm/sec  for  conventional  ignition.  The  equation  obtained  by  Zaghloul  et  al  [5]  yields  a 
bum  rate  of  33.72  cm/sec  with  plasma  injection.  In  general,  the  new  modified  form  follows  the 
same  trend  even  though  the  pressure  exponent  is  not  as  large.  This  is  due  to  the  large  difference 
in  the  magnitude  of  the  constant  ‘a’  and  the  pressure  range  looked  at. 
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The  equation  in  this  form  does  provide  a  good  application  oriented  equation  because  impact 
pressure  is  a  common  measurable  value  in  cannon  ballistics. 

It  is  also  possible  to  remove  the  pressure  dependence  by  taking  the  ratios  of  the  bum  rates  at  a 
raised  bed  temperature  to  the  bum  rate  at  ambient  bed  temperature.  Thus,  the  equation  may  then 
be  written  as: 


(BR\  aPbTr 
(BR\  ~  aPbTl 


=d 


f  T  \ 

\ToJ 


r 


JA2 


where  C  is  a  constant.  This  equation  is  valid  for  other  ignition  means,  e.g.  conventional  or  laser, 
because  it  depends  only  on  the  propellant  bed  temperature. 

Each  bum  rate  data  set  was  analyzed  and  was  normalized  to  the  ambient  propellant  bum  rate. 
The  values  used  for  the  normalization  were  the  bum  rates  at  T=20°C  or  the  average  of  the  coolest 
ambient  data  point  available.  The  1-inch  data  was  normalized  to  34.325,  2  inches  to  8.566,  and 
2.75  inches  to  3.485  cm/sec.  Normalized  bum  rate  is  shown  in  Fig.  13. 


Fig.  13  Normalized  bum  rate  of  JA-2  with  increased  bed  temperature  for  all  samples’  positions. 
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It  can  be  seen  from  Fig.  13  that  the  normalized  bum  rates  at  1  and  2  inches  are  almost  the 
same  while  the  normalized  bum  rate  at  2.75  inches  increases  significantly. 

The  normalized  equation  was  transformed  into  a  log  form  and  a  linear  fit  was  performed. 


In 


r(BR)T ' 


(BR) 


=  ln(C)  +  y  In 


To  J 


f  rp  \ 

V^O  Jjai 


The  fit  parameters  are  shown  in  Table  3  for  1,  2,  and  2.75  inches  from  the  electrothermal  source 
exit. 


Table  3:  Temperature  Exponent  y  and  constant  C  for  1, 2  and  2.75  inches 


C 

Y  . 

1-inch 

0.6110 

1.9721 

2-inches 

0.5638 

2.4767 

2.75-inches 

1.0352 

3.3312 

These  fits  are  shown  in  Figures  14,  15  and  16. 


Fig.  14  Normalized  Bum  Rate  Fit  for  d=l  inch  from  the  source  exit 
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Fig.  1 5  Normalized  Bum  Rate  Fit  for  d=2  inch  from  the  source  exit 


big.  16  Normalized  Bum  Rate  Fit  for  d— 2.75  inch  from  the  source  exit 


These  figures  show  that  the  propellant  bed  temperature  augments  the  data  more  the  farther  it 
extends  beyond  the  barrel,  which  indicates  the  bum  rate  is  also  dependent  on  the  distance  from  the 
source.  The  conclusion  here  is  that  the  plasma  parameters  do  play  an  important  role  in  propellant 
bum  and  need  to  be  incorporated  into  the  equations.  An  additional  functional  form,  with 
dependence  on  bed  temperature  and  distance  from  the  source  section  is  proposed. 
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BR 

BR, 


=  C 


'T\'  (dV 


\ToJjA2 


\doJ 


where  the  naught  values  are  taken  at  an  arbitrary  distance  from  the  source. 

Since  P  =  £  nkTp,  the  pressure  may  then  be  substituted  in  the  bum  rate  equation  to  include 
plasma  parameters  in  lieu  of  the  pressure: 


BR  =  a(nkT  )' 


r  T  v 


T 

V  oJjA2 


where  n  is  the  plasma  density,  k  is  the  Boltzmann  constant,  and  Tp  is  the  plasma  temperature. 

Plasma  temperature  can  be  related  to  plasma  heat  flux  utilizing  the  Blackbody  assumption 
such  that  the  plasma  temperature  could  be  expressed  in  terms  of  the  blackbody  heat  flux  q" : 

1 


where  s  is  the  emissivity  and  a  is  the  Stefan-Boltzmann  constant.  For  simplicity,  the  emissivity  s  is 
assumed  to  be  equal  to  1,  and  thus,  substituting  for  the  plasma  temperature  yields  an  expression 
for  the  bum  rate  in  terms  of  radiation  heat  flux: 


BR  =  Anb 


T 


V 


T 


O  J 


where  A  is  a  constant. 

Optical  emission  spectroscopy  has  been  performed  to  determine  the  plasma  temperature  and 
density.  Copper  lines  were  used  to  obtain  plasma  temperature  and  number  density.  The  lines  used 
for  the  temperature  analysis  and  their  parameters  are  shown  in  Table  4. 
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Table  4:  Copper  lines  used  for  Spectroscopic  Analysis 


Spectrum 

Wavelength  (nm) 

E,  (eV) 

gA(10V) 

Cu  I 

510.554 

3.82 

0.080 

Cu  I 

515.324 

6.19 

2.4 

Cu  I 

521.820 

6.19 

4.5 

Cu  1 

529.250 

7.74 

0.436 

Cu  I 

570.020 

3.82 

0.0096 

Cu  I 

578.213 

3.79 

0.033 

Three  spectral  samples  are  shown  for  each  distance  in  Figures  17,  18  and  19,  to 
demonstrate  reproducibility  of  the  spectral  data,  thus  indicating  similar  properties. 


Fig.  17  Sample  three  d  =  1-inch  shot  Spectra 
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The  PEAKFIT  [29]  program  was  used  to  determine  all  measured  peaks  in  the  spectra.  It  can 
determine  all  the  peaks,  obvious  and  hidden,  and  outputs  them  and  their  areas  to  file  for  analysis. 
Table  5  shows  plasma  temperature,  density  and  the  corresponding  radiation  heat  flux  at  the  three 
locations  1,  2  and  2.75  inches  from  the  source  exit. 


Tab 

e5:  Average  Plasma  Temperature,  Heat  Flux, 

and  Density 

Distance 

Temperature 

Heat  Flux  *  1 08 

Density  *1024 

_ (in) 

(eV) 

W/m2 

-3 

m 

1 

0.846 

5.237 

3.575 

2 

0.687 

2.289 

3.406 

2.75 

0.821 

4.666 

2.614 

The  purpose  for  determining  these  temperatures  is  to  attempt  to  determine  the  bum  rate 
dependence  on  plasma  parameters  as  well  as  propellant  bed  temperature.  A  bum  rate  ratio  can  be 
determined  between  the  distances  that  will  allow  correlation  with  number  density  and  plasma  heat 
flux.  The  ratio  can  be  expressed  as: 


(BR) 
( BR)a 


=  A 


\no  J 


(  v  Y 


f  j1  V 


lTy,,2 


The  only  question  about  this  equation  is  what  the  normalizing  condition  should  be.  The 
reasonable  choice  is  the  plasma  density  because  it  dominates  the  plasma  temperature  in  the  kinetic 
pressure  term  P  =  £  nkT.  For  this  reason,  all  naught  values  are  those  of  the  2.75  inch  data  at 
ambient  temperature  so  that  the  bum  rate  ratio  would  reduce  to  one  at  that  point.  Utilizing  the 
values  obtained  for  plasma  temperature  and  density  allows  the  unknowns  of  the  above  equation  to 
be  determined  using  Matlab’s  nonlinsq  function  and  are  shown  in  Table  6. 


Table  6:  Parameters  of  the  modified  equation  with  dependence  on  plasma  radiation  heat  flux  and 
_ number  density,  and  the  normalized  bed  temperature  with  respect  to  ambient _ 


Parameter 

- .  A _ _ _ 

Value 

A 

1.8315 

b 

2.6114 

_ 1 _ 

2.8365 

Hence,  the  bum  rate  equation  can  now  be  written  as: 
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This  Equation  provides  an  insight  into  the  effect  of  plasma  density  and  radiation  heat  flux  on 
the  bum  rate  of  the  propellant  bed.  It  is  now  possible  to  predict  bum  rate  augmentations  based  on 
the  plasma  density  and  heat  flux  variations  as  well  as  propellant  bed  temperature. 


IV.  Conclusion 

The  plasma  flow-field  parameters  were  measured  showing  the  distribution  of  essential 
parameters,  such  as  pressure,  temperature,  number  density  and  average  speed.  The  importance  of 
such  plasma-flow-field  characterization  is  to  help  in  the  design  of  electrothermal-chemical  devices 
by  knowing  the  distribution  of  the  plasma  jet  parameters. 

It  is  concluded  that  there  is  a  power  law  bum  rate  dependence  on  propellant  bed’s 
temperature.  This  has  been  shown  experimentally  and  has  a  reasonable  physical  meaning  because 
of  the  extra  energy  given  to  the  propellant  bed  due  to  the  heat  energy  input.  The  functional  form 
of  bum  rate  equation  has  been  expressed  in  various  forms  but  always  includes  the  bed’s 
temperature.  A  useful  form  is  that  of  the  normalized  bed  temperature  with  respect  to  ambient 
temperature  and  a  dependence  on  the  pressure.  Another  form  is  the  ratio  form  between  bum  rates 
at  elevated  temperature  to  that  at  ambient.  The  form  that  includes  plasma  parameters  in  terms 
either  of  plasma  temperature  or  radiation  heat  flux  is  of  more  importance  in  order  to  characterize 
effect  of  plasma  augmentation  mechanism  on  bum  rates. 
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